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We present the results of a study of the charmless vector-vector decay 5° — > toK*'^ with 657 x 10^ 
BB pairs collected with the Belle detector at the KEKB e+e~ collider. We measure the branching 
fraction to be B{B" -> cjK*") = [1.8 ± 0.7(stat) ± 0.3(syst)] x 10"" with 3.0a significance. We also 
perform a helicity analysis of the lj and A'*° vector mesons, and obtain the longitudinal polarization 
fraction fhiB'^ — > u)K*°) = 0.56 ± 0.29(stat)_Qgg(syst). Finally, we measure a large non-resonant 



branching fraction B[B^ -^ ujK^-k 
10~® with a significance of 9.5(t. 



■,Mk7v e (0.755,1.250) GeV/c^] = [5.1 ± 0.7(stat) ±0.7(syst)] 



PACS numbers: 13.25.Hw, 12.15.Hh, 13.88.+e 



The study of branching fractions and angular distri- 
butions of B meson decays to hadronic final states tests 
our understanding of both weak and strong interactions. 
Recently, B decays mediated by 6 — > sqq penguin am- 
plitudes have received much attention in the literature. 
Unlike b —^ c spectator amplitudes (which are much bet- 
ter measured), penguin amplitudes contain an internal 
loop and thus are potentially sensitive to new propaga- 
tors and couplings. Such mediating particles may have 
an energy scale too high to access directly. Several mea- 
sured b -^ sqq decays have yielded unexpected results; 
e.g., the decays B -^ <t)K* and B -^ pK*^ are found to 
have large transverse polarization []|, and B decays to 
the closely related final states K^n~ and if +7r° exhibit 
different patterns of direct CP violation [4]. These re- 
sults are difficult to accommodate within the Standard 
Model and may indicate the presence of new physics [3| . 
Furthermore, b — > sqq decays are useful for determining 
the angles (/)2 and cj)^ of the unitarity triangle [Jj . 

In this Letter we present a study of the b — *■ sdd decay 
B^ -^ iuK*^. Theoretical calculations for the branching 
fraction cover the range (0.3 — 10.0) x 10"^ [5J. Pre- 
viously, this mode has been searched for by CLEO [6[ 
and BaBar [7|; the latter group observed an excess of 
events with a significance of 2Aa. Our analysis uses 
605 fb"^ of data containing 657 x 10^ BB pairs; this 
sample is almost three times larger than that used in 
Ref. [7|. With this large data set we are able to measure 
both the branching fraction and longitudinal polariza- 
tion fraction for B" -^ loK*^ , and the branching fraction 
for non-resonant _B" -^ loK'^'k~ . The data were col- 
lected with the Belle detector Qj at the KEKB e+e" 



asymmetric-energy (3.5 GeV on 8.0 GcV) collider with a 
center-of-mass (CM) energy at the T(4S') resonance. The 
production rates of B^B*^ and B^B^ pairs are assumed 
to be equal. 

The Belle detector is a largc-solid-angle spectrometer. 
It includes a silicon vertex detector (SVD), a 50-layer 
central drift chamber (CDC), an array of aerogel thresh- 
old Cherenkov counters (ACC), time-of-flight scintilla- 
tion counters (TOF), and an electromagnetic calorime- 
ter (ECL) comprised of CsI(Tl) crystals located inside a 
superconducting solenoid coil that provides a 1.5 T mag- 
netic field. 

The S-daughter candidates are reconstructed through 
the decays uj -^ 7t~^7t~7t'^, K*'^ — > K'^7t~ and 7r° -^ 
77 [lO|. A charged track is identified as a pion or kaon 
by using particle identification (PID) information from 
the CDC, ACC and TOF systems. We reduce the num- 
ber of poor quality tracks by requiring that \dz\ < 4.0 cm 
and dr < 0.2 cm, where \dz\ and dr arc the distances of 
closest approach of a track to the interaction point along 
the z-axis (opposite the direction of the positron beam) 
and in the transverse plane, respectively. In addition, 
we require that each charged track have a transverse mo- 
mentum pt > 0.1 GeV/c and a minimum number of SVD 
hits. Tracks matched with clusters in the ECL that are 
consistent with an electron hypothesis are rejected. 

Photons used for tt" reconstruction are required to have 
energies in the laboratory frame greater than 50 (100) 
MeV for the ECL barrel (endcap), which subtends 32° — 
129° (17° - 32° and 129° - 150°) with respect to the 
beam axis. We require tt*' candidates to have an invariant 
mass in the range M^^ € (117.8, 150.2) MeV/c^ (±3cr in 



rriT^a resolution) and a momentum in the laboratory frame 
p^o > 0.39 GeV/c. 

We select w mesons with an invariant mass in the 
range M^^^ e (0.730,0.830) GeV/c^, and K*° mesons 
with Mk^ e (0.755,1.250) GeV/c^. These windows in- 
clude sideband regions to provide discrimination between 
signal and background components in the maximum- 
likelihood (ML) fit described below. To reduce combi- 
natorial background arising from low-momentum kaons 
and pions, we require that cosS^*" > —0.8, where Oj^*o 
is the K*^ helicity angle defined as the angle between the 
direction of the K^ and the direction opposite to the _B" 
momentum in the K*'^ rest frame. 

Signal decays are identified using the energy differ- 
ence {AE) and the beam-energy-constrained mass (Mbc). 
These are defined as AE = Eb — -Ebeam and Mbc = 
\/^beam - Pb^ where i^beam dcuotes the beam energy 
and Eb and pb denote the energy and momentum, 
respectively, of the candidate B-meson, all evaluated 
in the e+e~ CM frame. We retain events satisfying 
\AE\ < 0.2 GeV and Afbc e (5.20,5.29) GeV/c^, and 
define a signal region AE e (-0.10,0.06) GeV, A/bc G 
(5.27,5.29) GeV/c2. 

The dominant source of background is continuum 
e'^e" —fqq{q=: w, d, s, c) production. To discriminate 
relatively spherical BB events from jet-like qq events, we 
use 16 modified Fox- Wolfram moments (combined into a 
Fisher discriminant J^ [111), the CM polar angle between 
the B direction and the z-axis (Ob), and the displacement 
along the z-axis between the signal B vertex and that of 
the other B in the event (AZ). Further discrimination is 
provided by a 6-flavor tagging algorithm [l3| , which iden- 
tifies the flavor of the B meson accompanying the signal 
candidate via its decay products: charged leptons, kaons, 
and A's. This algorithm yields a quality factor r, which 
ranges from zero for no flavor discrimination to one for 
unambiguous flavor assignment. 

We use Monte Carlo (MC) simulated signal [Hi 
and data sideband events [defined as Mbc G 
(5.20,5.26) GeV/c^, \AE\ < 0.2 GeV ] to obtain prob- 
ability density functions (PDFs) for JF, cos 9b and AZ. 
These are multiplied together to form signal (Cg) and 
qq background (C^) likelihood functions, and we require 
that TZqq = Cg/{Cg + C^-q) be above a threshold. We 
divide the events into six bins of r and determine the 
optimum TZqq threshold for each bin by maximizing a 
figure-of- merit S/y/S + B, where S (B) is the number 
of signal (background) events in the signal region. This 
optimization rejects 99% of the qq background while pre- 
serving 50% of the signal. 

The fraction of events having multiple candidates is 
12%. We choose the candidate in an event to be the one 
that minimizes the quantity jM^-,,- m^o |. From MC stud- 
ies we find that this choice selects the correct candidate 
90% of the time. We also find that 9.6% of signal decays 
have at least one particle incorrectly identified but pass 



all selection criteria; these are referred to as "self-cross- 
feed" (SCF) events. 

We obtain the yields using a four-dimensional (4D) 
extended unbinned ML fit to AE, Afbc, Mt^-^tt and Mk-k- 
The likelihood function is given by 
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where Yj is the yield of component j, P* is the PDF 
for component j, and i runs over all events in the sam- 
ple. We include PDFs for the signal, qq background 
(gg), charm i3-decay background (6 -^ c), charmless 
B-decay background (6 -^ s,u,d), and non-resonant 
B^ -^ ujK^tt^ decays. The MC acceptances for non- 
resonant B° -^ K*^TT+TT-TT° and B° -^ K+ tt' tt+ 11' tt° 
are negligibly small and thus we do not consider these 
channels. 

The PDF for each component is defined as V^j ~ 
V, [AE^JV, [MIJV, (M;^ JP, (M], J. For the signal and 
LuK^TT components, we split the PDFs into two parts: 

^j = (1 - /scF)'Ptruc + /scf'Pscf' ^here /scf is the 
SCF fraction (17% for u!K~^Tr~), and "true" represents 
the correctly reconstructed decays. For the qq, b ^ c 
and b — > s,u,d backgrounds, no sizable correlations are 
found among the fitted variables. For the signal and 
ujK'^t:'' components, there are small correlations that 
are accounted for as described below. 

The K*^ and u; resonances are modeled with Breit- 
Wigner functions whose widths are fixed to their 
PDG [IJ] values. The Breit-Wigner function used to de- 
scribe the to resonance is convolved with a Gaussian of 
a = 5.7 MeV to take into account the detector resolution. 
This value, along with the means for both resonances 
and the fraction of qq background events containing w's 
and fir*°'s, are obtained from fitting the Mk-k and Mtttttt 
spectra of events in the data sideband. 

All other PDF shapes are obtained from MC simu- 
lation. For the signal and loK^-k^ PDFs, the sum of 
a Crystal Ball line shape [l5| and Gaussian is used to 
describe AE, and the sum of two Gaussians is used to 
describe Mbc. To take into account small differences be- 
tween the MC simulations and data, the Mbc and AE 
shapes for the signal and wiiT+Tr" PDFs are corrected 
according to calibration factors determined from a large 
B^ -^ D^ p^ , D^ -^ K~^Tr~Tr~ control sample. The 
Mktt PDF for wA'+TT" decays is represented by a thresh- 
old function with parameters determined from MC events 
where the Kit final state is distributed uniformly over 
phase space. 

For the qq background, we use a threshold ARGUS |16j 
function to describe Mbc, and linear functions to de- 
scribe AE and the combinatorial shapes of Mtttttt and 
Mk-k- The Mbc and AE shapes of the b ^ c background 
are described by an ARGUS function and a second-order 
Chebyshev polynomial, respectively. The remaining PDF 
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FIG. 1: Projections of the fit results onto (a) AE, (b) Mbc, 
(c) Mtttt-w and (d) Mktt for candidates satisfying (except for 
the variable plotted) the criteria AS G (-0.10,0.06) GeV, 
Mbc G (5.27,5.29) GeV/c^ and Mk^ G (0.755, 1.050) GeV/c^ 
The curves are for ujK*^ (dashed), ojA'^tt" (dot-dashed), the 
sum of the backgrounds (dotted), and the total (solid). 



shapes are modeled with non-parametric PDFs using 
Kernel Estimation 



m 



The following parameters vary in our final fit to the 



data: the signal, ujK^ 



c and qq yields, and the 



qq PDF parameters describing the AE, Mbc and com- 
binatorial shapes of M-^TTTT and Mk-k- The fraction of 
h ^> s,u,d events {fb—>s,u.d) is small (1.6%) and fixed to 
the MC value. The /scf for signal and ujK^tt^ decays 
are also fixed to their MC values. 

The fit results arc listed in Tableland the projections 
are shown in Fig. [TJ With the fitted yields Y , we cal- 
culate the branching fraction B as Y/{eMC ■ epid ■ N^g), 
where £mc is the event selection efficiency including 
daughter branching fractions obtained from MC simu- 
lation, N^g is the number of BB pairs produced, and 
e-pYD is an efficiency correction for the charged track se- 
lection that takes into account small differences between 
MC values and data. Our signal MC simulation is gen- 
erated with /l = 0.5; the change in acceptance for other 
values of /l is taken as a systematic error. For wif+7r~, 
eMC and B are for Mk-k € (0.755,1.250) GeV/c^ The 
significance is defined as ^— 21n(£o/'Cmax), where £max 
(£o) is the value of the likelihood function when the yield 
is allowed to vary (set to 0). The systematic uncertainty 
is included by convolving the likelihood function with a 
Gaussian whose width is equal to the systematic error. 
For signal and ujK'^tt~ decays, we account for small cor- 
relations between the fitted variables by fitting ensembles 
of simulated experiments containing all signal and back- 
ground components. The correlations give rise to biases 



TABLE I: Signal yield Y and its statistical uncertainty, 
MC efficiency emc, PID efficiency £pid, significance S with 
systematic uncertainties included, and measured branching 
fraction B. For loK'^-k~ , emc and B are for Mkt: G 

„2 



(0.755, 1.250) GeV/c^ 
tical (systematic) 



For B, the first (second) error is statis- 



Mode 


Y 


EMC (%) 


epiD 


S 


B (10"") 


a;i^*" 


32.9+-- 


2.91 


0.94 


3.0 


1.8 ±0.7 ±0.3 


L0K^-K~ 


146.5+^i;-| 


4.66 


0.94 


9.5 


5.1 ±0.7 ±0.7 



of ±2.9 and ±8.1 events for signal and ujK^-k , respec- 
tively. We correct the fitted yields for these biases. 

The main sources of systematic error are: track recon- 
struction efficiency (1.2% per track); tfi efficiency (4%); 
PID (1.3%); TV^B (1.4%); MC statistics (0.6%); PDF 
shapes {t\l%)\ h^s,uA (2.6%); /scf (l5;2%); the AS 
fit range (lo:o%); fitting bias (7.8%); the effect of higher 
A'*" resonances (tao%); h {t\l%)\ and 7^,^ (2.8%). 
The errors on the PDF shapes are obtained by varying 
all fixed parameters by ilcr and taking the fractional 
change in the yield as the systematic error. To obtain the 
error due to /scf and fb^s,u,d, we vary these fractions 
by ±50%. The uncertainty in the yield bias correction is 
taken to be the sum in quadrature of the statistical uncer- 
tainty on the correction and half the correction value. We 
consider the effects of higher K*^ resonances by includ- 
ing a PDF for B° -^ ujK^{U30)° and repeating the 4D 
fit with the ujKg{lA30)'^ yield fixed to the value obtained 
by extrapolating from a higher Mktt region. The error 
due to the uncertainty in /^ is obtained by varying /l 
by its errors measured below. To obtain the uncertainty 
due to the 7?,gg requirement, we vary the TZqq thresholds, 
and we also calculate the data/MC efficiency ratio for 
the B^ -^ D^ p^ , D^ -^ K~^TT~Tr~ control sample. 

We study the effects of interference between i?° -^ 



LoK*'^ and i?° -^ ujK~^tt~ decays as follows. We modify 
the Breit-Wigner PDF describing the K*'^ resonance of 
the signal to include an interfering amplitude and phase; 
for lack of more information, we take this amplitude to be 




MK^(GeV/c1 

FIG. 2: Signal ± loK+tv~ yields obtained from 2D fits to AE 
and Mbc in bins of A4ktt- The curves are for ujK*° (dashed), 
u!K'^Tv~ (dot-dashed), and the total (solid). 



constant in Mktt- Wc uniformly vary the amplitude and 
phase from zero to a maximum and, for each case, gener- 
ate and fit a large ensemble of toy MC experiments. The 
rms spread of deviations about the true value is taken 
as the systematic error (tii% for wK*^). Combining 
all errors in quadrature gives a total systematic error of 
{-i4()%)- The systematic errors considered for ujK'^tt^ 
are similar; the total is (ti3.8%)- 

To verify the large luK^tt' contribution (see Table [l|, 
we bin the data in Mkt^ from 0.65 — 1.25 GeV/c^ and, 
for each bin, perform a two-dimensional (2D) fit to AE 
and Mbc. The likelihood function consists of three com- 
ponents: signal + LuK^TT^ , qq + b ^ c, and b -^ s,u,d. 
We plot the resulting yields of signal -I- uiK'^'k' as a 
function of Mktt (Fig- E]) and fit this distribution to ex- 
tract the signal and ujK'^tt~ components. For Mk-k € 
(0.755, 1.250) GeV/c2 we obtain yields of 29.3 ± 12.1 and 
161.5 ± 16.3 for signal and w-ftT+Tr", respectively; these 
values are in good agreement with the results of the 4D 
fit after accounting for the fit bias. 

The differential decay width, after integrating over 
the angle between the decay planes of the ui and K*^ 
mesons, is proportional to (1 — f L)siv? 9 ^^sm^ 9k 'O + 
AJlCOS^ 9^cos'^ 9k''0. Here, 9^ is the uj helicity angle 
defined as the angle between the normal to the three- 
pion decay plane and the negative of the i?° momentum 
in the lo rest frame. The fraction of longitudinal polar- 



ization fL = 1^0 r/E; 



I^aP, 



where A-,^ are the helicity 



amplitudes for the longitudinal (A = 0) and transverse 
(A = ±1) states [18|- To determine /l, we bin the data 
in I cos9t^\ and cos0/^*o and, for each bin, perform a 4D 
fit to AE, Mbc, Mktt, and Mtttttt. The resulting signal 
yields as a function of the helicity cosines are shown in 
Fig. [3l We perform a simultaneous x^ fit to these dis- 
tributions, where the only floating parameter is Jl- The 
PDFs for the Aq and A^-^ helicity states are determined 
from MC simulation to take into account the detection 
efficiency. The statistical error is obtained from a toy 
MC study (the rms spread of the residuals from a large 
ensemble), since the errors in the distributions of Fig. [3] 
are correlated. Using a large toy MC sample we mea- 
sure a 2% bias in the fitting procedure, which we use to 
correct the central value. 

There are six main sources of systematic error in fi^: 
uncertainty in the PDF shapes (-1-0.16,-0.06); the frac- 
tions fb^s,u.d (+0.02,-0.01) and /scf (+0.01,-0.01); 
fitting bias (+0.02,-0.00); interference (+0.02,-0.01); 
and the TZqq requirement (+0.08, —0.05). Adding the var- 
ious systematic contributions in quadrature, we obtain a 
longitudinal polarization fraction 



/i(BO -^ ujK*°) = 0.56 + 0.29(stat)+°;og(syst). 



(2) 



In summary, using 657 x 10^ BB pairs we have found 
evidence for the B^ -^ luK*'^ decay with a signifi- 
cance of 3.0(7. We measure the branching fraction to 
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FIG. 3: Signal yields obtained from 4D fits to AE, Afbc, Mk^, 
and A/tttttt, in bins of (a) | cos 6ui\ and (b) cosSj^.o. The solid 
histograms show the results of the simultaneous x^ fit- The 
dashed (dot-dashed) histograms are the Aq (^±i) component. 



be B{B° -^ ujK*°) ^ [1.8 ± 0.7(stat) ± 0.3(syst)] x 10"*^. 
Our result is in agreement with theoretical estimates [51] , 
and with the central value obtained by BaBar [7[. We 
also perform a helicity analysis of the uj and K*'^ vector 
mesons and measure a longitudinal polarization fraction 
/l(S° -^ ujK*°) = 0.56 ± 0.29(stat)+o;Jg(syst). This 
central value is lower than that predicted by most theo- 
retical models but is similar to that measured for other 
b — > sqq decays [l[. In addition, we measure a large non- 
resonant branching fraction B[B'^ -^ loK^'k~;Mktt G 
(0.755, 1.250) GeV/c2] = [5.1±0.7(stat)±0.7(syst)]xl0-f^ 
with a significance of 9.5cr. Assuming a uniform phase 
space distribution, this implies a branching fraction of 
[79t}J(stat) ± ll(syst)] x lO"*" over the whole region. 
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